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The density matrix expansion (DME) of Negele and Vautherin is a convenient tool to map finite- 
range physics associated with vacuum two- and three- nucleon interactions into the form of a Skyme- 
Uke energy density functional (EDF) with density-dependent couplings. In this work, we apply 
the improved formulation of the DME proposed recently in I arXiv :0910.4979l by Gebremariam et 
al. to the non-local Fock energy obtained from chiral effective field theory (EFT) two-nucleon 
(NN) interactions at next-to-next-to-leading-order (N^LO). The structure of the chiral interactions 
is such that each coupling in the DME Fock functional can be decomposed into a cutoff-dependent 
coupling constant arising from zero-range contact interactions and a cutoff-independent coupling 
function of the density arising from the universal long-range pion exchanges. This motivates a new 
microscopically-guided Skyrme phenomenology where the density-dependent couplings associated 
with the underlying pion-exchange interactions are added to standard empirical Skyrme function- 
als, and the density-independent Skyrme parameters subsequently refit to data. A Mathematica 
notebook containing the novel density-dependent couplings is provided. 

PACS numbers: 21.10.Dr, 21.60.Jz, 21.30.cb, 71.15.Mb 

Keywords: Density matrix expansion, non-empirical energy density functional 



> 
O 

in 

rn 
O 

o 



I. INTRODUCTION 



A longstanding challenge of nuclear theory is to cal- 
culate properties of nuclei starting from the vacuum 
two- and three-nucleon interactions. While impressive 
progress has been made in extending the limits of ab- 
initio methods beyond the lightest nuclei [iHll, the nu- 
clear energy density functional (EDF) approach remains 
the most computationally feasible method for a compre- 
hensive description of medium and heavy nuclei . Mod- 
ern parameterizations of empirical Skyrme and Gogny 
EDFs provide a good description of bulk properties and, 
to a lesser extent, of certain spectroscopic features of 
known nuclei. However, the lack a solid microscopic foun- 
dation often leads to parameterization-dependent pre- 
dictions away from known data and makes it difficult 
to develop systematic improvements. Fueled by interest 
in the coming generation of radioactive ion beam facili- 
ties, along with studies of astrophysical systems such as 
neutron stars and supernovae that require controlled ex- 
trapolations of nuclear properties in isospin, density, and 
temperature, there is a large effort currently underway 
to develop energy functionals with substantially reduced 
errors and improved predictive power, e.g. see Ref. 

One path forward focuses on empirically improving the 
analytical forms and fittingprocedures of existing phe- 
nomenological functionals |6l4lll|. In the present work. 
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we pursue a complementary approach that relies less on 
fitting empirical functionals to known data, but rather at- 
tempts to constrain the analytical form of the functional 
and the values of its couplings from many-body pertur- 
bation theory (MBPT) and the underlying NN and NNN 
interactions (13 - ll8l |. 

Recent progress in evolving chiral effective field theory 
(EFT) interactions to lower momentum using rcnormal- 
ization group (RG) methods (see also [H ^) 

plays a significant role in this effort, as the many-body 
problem formulated in terms of low-momentum interac- 
tions is simplified in several key respects. The evolu- 
tion to low-momentum weakens or largely eliminates non- 
perturbative behavior in the two-nucleon sector arising 
from strong short-range repulsion and tensor forces from 
iterated pion exchanges [H, [2^. In addition, at lower 
cutoffs the corresponding three-nucleon interactions be- 
come perturbativc and more amenable to approximations 
such as truncations based on normal-ordering pTl [28| . 
When applied to nuclear matter, many-body perturba- 
tion theory for the energy appears convergent (at least 
in the particle-particle channel), with calculations that 
include all of the NN and most of the NNN second-order 
contributions, exhibiting reasonable saturation proper- 
ties and showing relatively weak dependence on the cut- 
off [l^ll^,!!^. Moreover, the freedom to vary the order 
of the input EFT interaction and the cutoff via the RG 
provides a powerful tool to assess theoretical errors aris- 
ing from truncations in the Hamiltonian and many-body 
approximations. 

All of these features are favorable ingredients for the 
microscopic construction of non-empirical EDFs [Tsj . In- 
deed, Hartree-Fock becomes a reasonable (if not quanti- 
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tative) starting point, which suggests that the theoreti- 
cal developments and phenomenological successes of EDF 
methods for Coulomb systems may be applicable to the 
nuclear case for low-momentum interactions. However, 
even with these simplifications, perturbative contribu- 
tions to the energy involve density matrices and prop- 
agators folded with finite-range interaction vertices, and 
are therefore highly non-local in both space and time. In 
order to make such functionals numerically tractable in 
heavy open-shell nuclei, it is desirable to develop sim- 
plified approximations expressed in terms of the local 
densities and currents. At lowest order in MBPT (i.e., 
Hartree-Fock) , the density matrix expansion (DME) of 
Negele and Vautherin [30| can be unambiguously ap- 
plied to approximate the spatially non-local Fock expres- 
sion^ as a generalized Skyrme functional with density- 
dependent couplings calculated from vacuum interac- 
tions. In the present work, we do so on the basis of 
a chiral NN interaction at next-to-next-to-leading-order 
(N^LO), while the extension to a chiral NNN interaction 
at the same order will be discussed in a separate paper. 

The non-trivial density dependence of the DME cou- 
plings is a consequence of the finite-range of the under- 
lying NN interaction and is controlled by the longest- 
ranged components. Consequently, the DME offers a 
path to incorporate physics associated with long-range 
one- and two-pion exchange interactions into existing 
Skyrme functionals. Given the rich spin and isospin 
structure of such interactions, it is hoped that their in- 
clusion will improve predictive power away from known 
data and provide microscopic constraints on the isovector 
structure of nuclear EDFs. 

Still, calculations of infinite nuclear matter (INM) [22| 
as well as binding energies and charge radii of doubly- 
magic nuclei [l^l demonstrate that it is necessary to go 
at least to second-order in perturbation theory to re- 
sum enough bulk correlations. Furthermore, it is known 
that while chiral EFT interactions are themselves low- 
momentum interactions relative to conventional force 
models, it is still desirable to evolve them to lower- 
momentum so that HE becomes a reasonable starting 
point and MBPT is under better control. 

However, in the present paper we focus on the lowest- 
order (i.e., Hartrec-Fock) contribution to the energy from 
the un-evolved chiiail EFT NN interaction. In light of the 
proceeding remarks, this might appear to be an unrealis- 
tic starting point. This would certainly be the case if our 
present goal was to develop a fully microscopic and quan- 
titative EDF free from any fitting to data. In the short 
term, however, we adopt a more pragmatic approach. 
Our objective in the present approach is to improve ex- 



isting Skyrme phenomenology by identifying non-trivial 
density dependencies arising from missing pion physics 
that can be added to existing Skyrme functionals, which 
can then be refit to data and implemented in existing 
codes with minimal modification. The rationale for re- 
stricting our attention to the Hartree-Fock energy using 
un-evolved chiral NN interactions can be summarized as 
follows: 

• First, it is well-known that the RG evolution to low 
momentum only modifies the short-distance struc- 
ture of the inter- nucleon interactions 0, [2^, [3l| . 
The input chiral NN interactions take the schematic 
form 



(1) 



^ We assume local NN interactions since our focus here is on the 
finite-range pion exchanges, which are local up to an overall cut- 
off regulator. For non-local interactions, the Hartree contribu- 
tion is no longer local in space in the sense that it probes the 
off-diagonal part of the density matrix. 



where denotes the finite-range pion-exchange 
interactions and Vj,^^(A) denotes scale-dependent 
zero-range contact terms. The RG evolution only 
modifies V^^ and leaves the long-distance struc- 
ture unchanged. Since we arc primarily interested 
in identifying the dominant density dependencies 
arising from finite-range physics, it is sufficient for 
our purposes to apply the DME directly to the un- 
evolved V^^ from the input EFT. 

Note that the HE energy arising from V™ bears a 
strong resemblance to the empirical Skyrme func- 
tional (i.e., bilinear products of local densities mul- 
tiplied by coupling constants), and therefore docs 
not produce any new density dependencies. 

• Second, a non-trivial extension of the DME is 
needed to treat non-localities in both space and 
time that arise in higher orders of perturbation 
theory. I.e., one must properly account for the 
presence of energy denominators when designing a 
DME for 2nd-order MBPT and beyond [H. To 
date, a satisfactory generalization of the DME has 
not yet been formulated. 

• Third, even if we were to follow the ad-hoc prescrip- 
tion that consists of replacing the vacuum NN inter- 
action in the Hartree-Fock expression by a Brueck- 
ner G-matrix (or a perturbative approximation in 
the case of low-momentum interactions) evaluated 
at some average energy, the G-matrix differs from 
the NN potential only at short distances. 

Therefore, while a Hartree-Fock calculation using the 
un-evolved chiral EFT NN interaction would provide a 
very poor description of nuclei, the application of the 
DME to such contributions captures some of the same 
density dependencies that would arise from the finite- 
range tail of any in- medium vertex (e.g., a G matrix 
or a perturbative approximation thereof) that sums lad- 
der diagrams in a more sophisticated many-body treat- 
ment. Once a satisfactory generalization of the DME is 
developed to handle spatial and temporal non-locality on 
the same footing, non-localities arising from in-medium 
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EDF Energy density functional 

DME Density matrix expansion 

PSA Phase space averaging 

NV Negele and Vautherin 

OBDM One-body density matrix 

INM Infinite nuclear matter 

MBPT Many-body perturbation theory 

HP Hartree-Fock 

EFT Effective field theory 



TABLE I: List of acronyms repeatedly used in the text. 



propagation can be mapped into the density-dependent 
Skyrme couplings as well^. 

The rest of the paper is organized as follows. In Sec.HIl 
we derive the Hartree-Fock energy for even-even nuclei, 
which serves as the starting point for the DME. Sec- 
tion imi reviews the improved PSA-DME of Ref. [H 
that is used in the present work. Master formulas and 
skeleton expressions for the resulting DME couplings ob- 
tained from the Fock energy are given in Sec. lIVI Results 
for various density-dependent couplings arc discussed in 
Sec. |V] and conclusions are given in Sec. IVIl Various 
technical details and lengthy expressions are given in the 
appendices. The explicit forms for the chiral EFT NN 
finite range and contact interactions are given in Appen- 
dices |X] and |B1 and the exchange interaction is given in 
Appendix O The PSA-DME is reviewed in Appendix |D] 
Formulas to construct the single particle fields obtained 
from the density-dependent DME couplings are given in 
Appendix [El and the couplings that arise from perform- 
ing the DME on the Hartree energy are collected in Ap- 
pendix |F1 Finally, detailed expressions of the DME cou- 
plings are provided in a companion Mathematica note- 
book, and are also shown in Appendix [G] 



II. HF ENERGY FOR EVEN-EVEN NUCLEI 

Before applying the DME to the Fock contribution 
from the chiral NN interaction, it is useful to provide a 
detailed expression of the Hartree-Fock potential energy, 
Vhf- We restrict our attention to the ground states of 
even-even nuclei throughout this paper, with the con- 
sequence that certain contributions to the energy are 
zero due to the intrinsic time-reversal invariance of such 
states. Note however that the PSA-DME of Ref. [H 
provides a natural framework to extend the approach to 



^ It remains to be seen if such a generalization will necessitate the 
introduction of orbital-dependent terms into the EDF. 



states that break time-reversal symmetry, i.e. ground 
states of odd-even or odd-odd nuclei, see Appendix |D] 
For a general (possibly non-local) two-nucleon potential 
yNN ^ T/jjp is defined in terms of occupied self-consistent 
HF orbitals as 

1 ^ 

^HF = ^Y.{^j\V''''{l-Pr2m 
ij 

^ ^E(uIV^^Nj). (2) 

ij 

The antisymmetrizcd interaction V = 1/^^(1 — P12) has 
been introduced, with P12 equal to the product of spin, 
isospin and space two-body exchange operators P12 = 
PcrPrPr, whcre 

= i(l + ai.(72) and = ^(1 + n -tz) . (3) 

By making repeated use of the completeness relation 

1 = /dr |rgr)(ro-r| , (4) 

and the definition of the HF density matrix 
A 

p{r3<73T3,ri(TiTi) = E(/i*(ri(TiTi)0,(r3cr3T3) , (5) 

i 

Eq. ([2]) can be written as 
1 f * 

Vhf = 2 X! / n^''* ^'■i'^i'^i''2cr2-r2|V^^|r3cr3r3r4(T4r4) 

{ary i=l 

X p(r3Cr3T3, riaiTi)p{r4a4,T4, V2a2Ti) 

= ^TriTr2/'ndr.(rir2|Vj^^2|r3r4) 

xp(i)(r3,ri)p(2)(i.4,r2), (6) 

where a matrix notation in spin and isospin spaces is 
used in the second equation and the traces denote sum- 
mations over spin and isospin indices for "particle 1" and 
"particle 2". Hereafter we drop the indices on "y^^ and 
p indicating which space they act in as it will be clear 
from the context. Switching to relative/center-of-mass 
(COM) coordinates and noting that the free-space two- 
nucleon potential is diagonal in the COM coordinate, the 
Hartree-Fock expression becomes 

Vhf = iTriTr2y"dRdrdr'(r'|V^^|r) 

x p(R + r/2,R + r72) 

xp(R-r/2,R-r72), (7) 
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where the antisymmetrized coordinate space interaction 
is given by the Fourier transform 

I 



where p' and p are the "incoming" and "outgoing" rela- 
tive momenta and (pIV^'^Ip') is understood to be an op- 
erator with respect to spin/isospin quantum numbers and 
a matrix element with respect to momentum. For chiral 
NN interactions through N^LO, the explicit spin/isospin 
structure can be expressed as 



(p'|V^^|p) = [Vc + Ti ■T2Wc] + [Vs + Ti ■T2Ws] 0-1 •a-2 + [Vt + Ti ■T2Wt] ai ■qo-2 -q 

+ \Vc + Tl • T2 Wc] + [ Vs + Tl ■ T2 Ws ] CTl ■ (72 + [Vt + Tl ■ T2 Wt ] CTl ■ k CT2 ' k 

+ [Vis + Tl ■T2 W^Ls] ^ (o-l + (72) • (q X k) , (9) 

where the two-body form factors {Vc, Vs, ■ • ■} are functions of the momentum transfer q = p' — p and {Vc, Vs, • • •} 
are functions of k = (p + p')/2. Explicit expressions for the momentum-space form factors are given in Appendix VK\ 
The k-depcndent terms arise entirely from the zero-range contact terms in the chiral EFT potential, while the 
q-dependent terms receive contributions from both finite-range pion exchanges and zero-range contact terms. De- 
composing (p'l V^^P'^'^I — p) analogously to Eq. [51 one has 

(p'|y^^P""| -p) = [vs + Ti-T2W^] + [V^ + n -TsM^l] (71 •a2 + [V^ + Ti-T2W^]ai-ka2-k 

+ + Tl • T2^V^] + [V^ + Tl ■ T2W^] ai ■ a2 + [vs + Tl • T2 VV|] (71 • q (72 • q 

+ [VSs + n-T2Wls]^{<Ji + <y2)-{kxq), (10) 

I 



where the exchange force form factors {VS, Vg, • ■ •} are 
functions of k whereas {VS, V^ , ■ • ■} are functions of q. 
The exchange form factors can be expressed as linear 
combinations of the direct interaction, and are given in 
Appendix 

Momentum space NN potentials that depend only on q 
or k are sometimes referred to as quasi-local, as their co- 
ordinate space representation is diagonal such that the 
anti-symmetrized interaction in Eq. [B] is schematically 
given by 

(r|V^"|r') = S{r - r')V{r, V) - 6{r + r')V{r, \7)P<'^ . 

(11) 

For the chiral potentials considered here, the V{r, V) can 
have terms that are linear (via the spin-orbit interaction) 
or quadratic (via the contact interactions) in V. In order 
to preserve the quasi-local nature of the potential, we ne- 
glect the ultra-violet (UV) regulator that multiplies each 
term in Eq. [Sj The most commonly used regulators give 
a non-local interaction since (p'|V^^^|p) in Eq. [H]is re- 
placed by f{p'/A){p'\V''''\p)fip/A), where f{p/A) ^ 
for p ^ A and f{p/A) « 1 for p ^ A. While it is possible 
to use a regulator that suppresses large momentum trans- 
fers instead of large relative momenta, it is reasonable to 
neglect the regulator since we work at the HE level with a 
local fermi momentum kp ^ A. Note that the UV cutoff 
A is between ~ 2.5 — 3.0 fm~^ in most implementations 
of the N^LO NN interactions. 



Given that we are working with a quasi-local interac- 
tion, it is convenient to treat Hartree and Fock contribu- 
tions to Eq. [7] separately. To do so, we first expand the 
p matrix on Pauli spin and isospin matrices 

P(ri,r2) = ^[po(ri,r2) +Pi(ri,r2)T;, 

+ So(ri,r2)-(7-hSi(ri,r2)-(7T^] , (12) 

where the usual scalar-isoscalar, scalar-isovector, vector- 
isoscalar, and vector-isovector components are obtained 
by taking the relevant traces, 

po(ri,r2) = Tr^^[/j(ri,r2)] , (13) 

pi(ri,r2) = Tiw[p(ri,r2)T,] , (14) 

So(ri,r2) = Tiw[p(ri,r2)a] , (15) 

Si(ri,r2) = Tiw[p(ri,r2)(7T,] . (16) 

Starting from these non-local densities, it is useful to 
define the following set of local densities (i = 0, 1) 



Pt(R) = Pt(ri,r2)|ri=r2=R , 
Tt(R) EE Vi • V2 /9t(ri,r2)|ri=r2=R 



(17) 
(18) 



Jt(R) = --(Vi-V2)x St(ri,r2)|ri=r.=R ,(19) 

which correspond to the matter density, the kinetic den- 
sity and the spin-orbit current density, respectively. 
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Inserting Eqs. |9] and [10] into Eq. [71 evaluating the 
spin/isospin traces and dropping terms that vanish in 



even-even nuclei (e.g., terms involving the local part of 
the spin density S(r,r)) finally gives 



and 



t=o,i 



= - ^ / dRrfrLt(R + r/2)pt(R-r/2)r*,(r)+r.Jt(R + r/2)pt(R-r/2)rl5(r) 



Vf = / f^R*[p?(R + r/2,R.-r/2)rg*(r) -S2(R + r/2,R-r/2)r;^*(r 

t=o,i •' 

+ Sf{K + r/2, R - r/2)5f (R + r/2, R - r/2) V„V0rf,*(r) 
+iY. I rfrirfr2rfs(r)St(r2,ri)- (r X Vi)pt(ri,r2), 

I 



t=o,i • 



(20) 



(21) 



where Vr and Vp denote the direct (Hartree) and ex- 
change (Fock) contributions, respectively. The F vertices, 
which in fact only depend on the norm of r, are defined 



F*(r) ^ y.(r)-lf(r) t = 0, 
= W,(r) - Wf(r) t = 1 , 

Ff(r) s Vr{v)^V^v) t = 0, 



(22) 



(23) 



for i S {C, 5, T, where the coordinate-space inter- 
actions are given by, e.g.. 



(27r)3 
I f dq 



e'^'V^iq) fori = C,S,T, (24) 



re^i"- (q-r)1^,(q) fori = LS.(25) 



As discussed in Rcf. [3J] , our primary focus is to apply 
the DME to the exchange (Fock) part of the HF energy 
while treating the Hartree term exactly. Indeed, it was 
realized long aga^tarting with the early works by Negele 
and Vautherin [sO, [131 , that 

(i) Treating the Hartree contribution exactly provides 
a better reproduction of the density fluctuations 
and the energy produced from an exact HF calcu- 
lation m. 

(ii) Restricting the DME to the exchange contribution 
significantly reduces the self-consistent propagation 
of errors [sj, [IHl • 

(iii) Treating the Hartree contribution exactly generates 
no additional complexity in the numerical solutions 
of the resulting self-consistent HF equations [131 
compared to applying the DME to both Hartree 
and Fock terms. 



Nevertheless, it is possible to apply the DME to the 
Hartree terms so that the complete Hartree-Fock con- 
tribution is mapped into a local EDF. For completeness, 
the DME couplings arising from the Hartree contribu- 
tions are also collected in Appendix [Fj However, owing 
to the deficiencies of the DME at reproducing such con- 
tributions, it is strongly advised to treat the finite range 
Hartree terms exactly in actual self-consistent EDF cal- 
culations 36 1 . See, however, Ref. [l^] where an accurate 
local EDF approximation for the Hartree energy for the 
Gogny force is obtained by performing a simple Taylor 
series expansion. 



III. DENSITY MATRIX EXPANSION 

Long ago, Negele and Vautherin formulated the den- 
sity matrix expansion to establish a theoretical connec- 
tion between the phenomenological Skyrme energy func- 
tional, which is written as sums over bilinear products 
of local densities and currents, and microscopic Hartree- 
Fock expressions involving the non-local density matrix 
and finite range NN interaction [s^ [33|. The central 
idea is to factorize the non-locality of the one-body den- 
sity matrix (OBDM) by expanding it into a finite sum of 
terms that are separable in relative and center of mass 
coordinates. Adopting notations similar to those intro- 
duced in Ref. |38[, one writes 



pt(ri,r2) « J2^n(kr)Vn{R) , (26) 

"I max 

St(ri,r2) « J2 ^^nikr) Q™(R) , (27) 



m=0 



where A: is a yet-to-be-specified momentum that sets 
the scale for the decay in the off-diagonal direction, 
whereas H{(fcr) denote the so-called H— functions that 
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also remain to be specified.^ Functions {■p„(R), Qm(R)} 
denote various local densities and their gradients 
{pt (R) , n (R) , Jt (R) , Vpt (R) , Apt (R) } . 



The benefit of expansion [26ll27l is to approximate the 
non-local Fock energy (Eq. [^T|) as a bilinear local Skyrme- 
like EDF of the form (for time-reversal invariant systems) 



J 



n 1 



t=o,i ' 



c-r^i-vpt 



(28) 



r 



where the densities depend locally on R, while the cou- 
plings that are microscopically derived from the vacuum 
NN interaction depend on the arbitrary momentum scale 
k. In the present work, we adopt the usual LDA choice 
where k is chosen to be the local Fermi momentum re- 
lated to the isoscalar density through 



fcF(R) 



Po(R) 



1/3 



(29) 



although other choices are possible that include addi- 
tional T- and Ap-dependencies [s^. The DME cou- 
plings are therefore density-dependent (or equivalently 
R-dependent) and are given by integrals of the finite- 
range NN interaction over various combinations of the 
H-functions, e.g. 



Cr(R) ^ / drT'IiP,{kFr)Ii^,{kpT)Tf{r), (30) 



where the R-dependence oi kp is suppressed for brevity. 
In Eq.[28l only pseudovector contributions to the more 
complete "tensor terms" [ioj have been kept for simplic- 
ity. While this is exact in spherical nuclei, it is only 
approximate in nuclei that break rotational invariance. 
Howev er, p seudotcnsor contributions have recently been 
shown [40] to be systematically two orders of magnitude 
smaller than vector ones in axially deformed nuclei, which 
justifies the common practice of neglecting the former for 
the purpose of calculating binding energies in situations 
where Galilean invariance is not broken. 



Several DME variants have been developed in the 
past [13, m, m. Ell, El]. They mainly differ regarding 
(i) the choice of the momentum scale fc, (ii) the path 
followed to obtain actual expressions of the 11— functions 
(see below) and (iii) the set of local densities that oc- 
cur in the expansion. As discussed in Ref. [s^l, all of 
these variants give reasonably accurate descriptions of 
the Fock energy contributions that probe the scalar part 
of the OBDM. 

However, for the (spin) vector part of the OBDM, 
which is relevant for approximating Fock contributions in 
spin-unsaturated nuclei where at least one pair of spin- 
orbit partners is only partially filled, the situation is very 
different. In Ref. j33| . it was shown that the poor accu- 
racy of the original DME of Negele and Vautherin (NV- 
DME) for the vector part of the OBDM can be dramat- 
ically improved by using phase space averaging (PSA) 
techniques and by accounting for the anisotropy of the 
local momentum distribution that is a generic feature of 
finite Fermi systems psj . The anisotropy is especially 
pronounced in the nuclear surface region where the vec- 
tor part of the OBDM is sharply peaked, and should be 
accounted for if the DME is to provide an accurate de- 
scription of spin-unsaturated nuclei that constitute the 
majority of nuclei. 

In this paper we present results for both the standard 
NV-DME [30| and a simplified variant of the recently 
developed PSA-DME [s^. For both versions, we follow 
common practice and truncate the DME to rimax = 2 in 
Eq. [2S] for the scalar part 



J 



,R 



^Apt(R)-Tt(R) + ^4pt(R) 



(31) 



and to TTlniax 

St(R- 



1 in Eq. [27] for the vector part 



,R 



-) - 
2' ~ 2' 2 

with the additional feature that the n = 1 (m 
contribution to the scalar (vector) part of the OBDM 
is zero in time-reversal invariant systems. The NV-DME 



-m{kFT)v X Jt(R), (32) 



0) 



H-functions arc given by. 



.7i(/cF(R)r) 



H^(fci.r) = 105 



/ci.(R)r ' 
j3(fcF(R)r) 
(fcF(R)O' 



^\{kFr) = jo(fcF(R)r), 



(33) 

(34) 
(35) 
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whereas the PSA-DME vector H-functions are given by 
ng(fcj^r) = nPikpr) = nUkFr) = 3 ^^^^^^^^^ ^ (36) 



kFi'R)r 



with 



(37) 



The function P2i'R) denotes the quadrupolc anisotropy 
of the local momentum distribution, which can be calcu- 
lated from the Husimi phase-space distribution'^ 



^ /rfp[3(e,-p)^-p^]g(r,p) 
/dpp2i7(r,p) 

A 

^|(e,-V)^,(r)p -1 



ro(r) ^ 



(38) 



where (pi (r) denotes an occupied HF single particle spinor 
with components ipiijcar). 

As noted in Ref. [s^ and discussed in Appendix |D1 
there is a simplified PSA approximation that uses the 
phase space of infinite nuclear matter and amounts to set- 
ting P2 (R) = in Eq. [5^1 From a practical perspective, 
the simplified PSA-DME provides a convenient compro- 
mise as it gives substantial improvements over the NV- 
DME in describing the vector part of the OBDM while 
avoiding the complicated single-particle fields that arise 
from the P2(R) dependence in the full PSA-DME. This 
is the version actually used in the present work. 



IV. DME COUPLINGS 



Separating long- and short-distance 
contributions 



Before evaluating the DME couplings, it is convenient 
to notice that a clean separation between long- and short- 
distance physics can be made due to the generic structure 
of the EFT interactions. Schematically, the EFT poten- 
tial can be written as 



V 



NN 



yr+T/,r+.-.+Kr(A) 



NN/ 



(39) 



where and are finite-range one- and two-pion 
exchange interactions dictated by the spontaneously bro- 
ken chiral symmetry of QCD, and V^^'^ denotes scale- 
dependent contact terms that encode the effects of inte- 
grated out degrees of freedom (heavier meson exchanges, 
high-momentum two-nucleon states, etc.) on low energy 



* The Husimi distribution _ff (r, p) is a positive-definite generaliza- 
tion of the Wigner /(r,p) function, see Ref. [43l . 



physics. Through N^LO, the contact interaction takes 
the form 

(Pl^rip') = + Ct ai • a2 + Ci q2 + C72 

+ "iCs ^ (o-i + CT2) • (q X k) + Ce (q • CTi)(q • CT2) 

+ C7(k-ai)(k-(T2) . (40) 

Note that the overall UV regulator has been neglected as 
discussed in Section|lIl and the A-dependence of the cou- 
plings is suppressed for brevity. For these contributions 
to the Hartree-Fock energy, there is no need to perform 
the DME since the zero-range nature of results in 
an expression that is already in the form of a bilinear 
Skyrme-like EDF with density-independent coupling con- 
stants (see Appendix IB]) . Consequently, each DME cou- 
pling at the HF level can be decomposed as the sum of 
a density-independent, A-dependent piece coming from 
Vj^^(A) and a density-dependent, A- independent piece 
coming from finite-range pion exchanges, i.e. 



cr 



Cr{A;VD + Cr{B.;Vr),etc. (41) 



As discussed in the Introduction, while Hartree-Fock 
becomes a reasonable zeroth-order approximation with 
low-momentum interactions, it is necessary to go to at 
least 2nd-order in MBPT to obtain a reasonable descrip- 
tion of bulk properties of infinite nuclear matter as well 
as binding energies and charge radii of closed-shell nu- 
clei. However, a consistent extension of the DME be- 
yond the Hartree-Fock level to treat the state-dependent 
energy denominators that arise has not been formulated 
to the best of our knowledge. At this point in time, 
any attempt to construct a quantitative Skyrme-like EDF 
starting from microscopic many-body theory must in- 
evitably resort to unsystematic approximations (e.g., av- 
erage state-independent energy denominators) in per- 
forming the DME on iterated contributions beyond the 
HF level [USHIilliil, and/or to the re-introduction 
of some phenomenological parameters to be adjusted to 
data. In the present paper, we are motivated by the fol- 
lowing observations 

1. A more quantitative many-body treatment such as 
2nd-order MBPT or the Brueckner-Hartree-Fock 
(BHF) approximation can be approximately cast 
into the same form as Eas.[^and[?T]if one neglects 
both the state-dependence of the intermediate- 
state energy denominators and the non-locality of 
the corresponding G matrix vertex. 

2. The in- medium G matrix vertex differs from the 
free-space NN potential mostly at short distances 
and "heals" to the free-space potential at suffi- 
ciently large distances. In the zeroth approxima- 
tion, this amounts to a fci?(R)-dependent renormal- 
ization of the couplings for the contact interaction 



NN 
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3. The coupling constants {Cg, Ct, Ci, . . .} of V^^ 
are usually matched to low-energy NN scattering 
data and deuteron properties, although in principle 
they could be matched to low-energy finite nuclei 
data. 



Based on these observations, we advocate a semi- 
phcnomenological approach in which the DME couplings 
from the finite-range Fock energy contributions arc added 
to empirical Skyrmc EDFs whose parameters are then 
re-fit to nuclear matter and finite nuclei properties. Al- 
though the treatment of the N^LO NNN contribution 
to the HF energy is postponed to a separate paper, the 
DME couplings from its long-range part constitute an 
integral part to be added to phcnomenological Skyrme 
parameters. Those re-fit parameters can thus be viewed 
as containing the effects of the HF contribution from the 
contact interactions V^^ plus higher order effects that 
would arise in a more sophisticated BHF or 2nd-order 
MBPT calculation. Finally, due to the loose connection 
of the refit Skyrme parameters to the EFT contact terms, 
the EFT concept of naturalness [l^] might provide use- 
ful theoretical constraints for the fitting procedure. The 
first calculations following this semi-phenomenological 
approach are underway and will be the subject of a future 
work 13611. 



IniQ/kp) functions are defined as 





f 2 

/ x'^dxjoiqx) (ng(x)) 


(46) 


h{q) ^ 


/ x'^dxjn{qx)n^{x)Uf^{x) 


(47) 


h{q) = 


j x^dxjoiqx) (n^(a;))^ 


(48) 


IM ^ 


1 xUxjoiqx) {nl{x)y. 


(49) 



Inserting the NV-DME and PSA-DME expressions for 
the H-functions gives 



^1(9) 
h{q) 
h{q) 

h{q) 



37r 
32 

35 TT 



t - 12q + 16) 61(2 - q) 



(50) 



128 



[f ~ 18q^ + 40g2 - 24g)6'(2 - q) (51) 



(52) 
(53) 



for the NV-DME and 



B. Master Formulas 

Inserting Eqs. l5T]and[5^into the Fock energy fEq. l^T]) 
and performing tedious but straightforward algebra to 
cast the expression into the same form as Eq. [28l we ob- 
tain the following "master formulas" for the DME cou- 
plings: 



cr 



'pAp _ 



a'-' = - 



1 



67rfc| 



h{q) = hiq) = — {f - 12q + 16) 9(2 - q) (54) 
h (q) = /4 (g) - ^ (2 - f) 0{2 - q) (55) 

oq 

for the PSA-DME, respectively. 

C. Skeleton Expressions 

The lengthy analytic expressions for the DME cou- 
plings obtained from the master formulas tend to ob- 
scure their underlying structural simplicity. Therefore, 
it is more illuminating to display the couplings in "skele- 
ton form" and relegate the explicit expressions to the 
Mathematica notebook and Appendix [Gl Each coupling 
Cf is given by the sum of the LO {n = 0), NLO (n = 1), 
and N^LO {n = 2) contributions 

2 

Cfiu) = 9 e P^P^ • ■ •} ' (56) 

n=0 

with u = /cf/wtt and 

2 

Cl^{u) =Y.'^'M^,u)J^^(n,u) , (57) 

where a^^(n,t,u) are rational polynomials in u and 
q''dqTg'^{q) l4{q/kF) , (45) Tj{n,u) are functions which may exhibit non-analytic 

behavior in u due to the finite-range of the NN interac- 
tion. 

We note that the detailed form of the skeleton NLO 
and N^LO expressions depends on the value of the spec- 
tral function regulator (SFR) mass used to regulate 



q'dqTfiq) [hiq/kp) 



+ ^hiqlkp)] (42) 



q^dqVf{q)h{q/kp) 



(43) 
(44) 



^ ^ q^dqh{q/kF){l+'^qd,)T-T\q) 



where terms with more than two gradients have been 
dropped, although investigating the effect of higher-order 
gradients might be of interest in the future [33|. The 
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the divergent loop integrals of the two-pion exchange po- 
tentials (TPEP), see Appendix H] and Ref. [H]. For sim- 
plicity, the following skeleton expressions were obtained 
for Mgfr oo, which corresponds to using dimensional 
regularization for the loop integrals that enter into the 
expressions for the NLO and N^LO two-pion exchange 
potentials as in Ref. [i^. In the skeleton expressions 
listed below, we use a more compact notation where the 
dependence on u, t, and n is not explicitly shown for the 
a's: 



LO 



= -f af log(l -I- iu^) + al arctan(2u) (58) 



NLO 



log(l + + 2u^/l + v?) 



+ al^/l + ^\og{l + 2v? + 2u\/T+v^) (59) 



a\ log(l + u^) + arctan(M) (60) 

Note that rather similar results are obtained for both the 
PSA-DME and NV-DME, as wcU as for finite values of 
the SFR mass Mgb in the TPEP. In the following section, 
we present results for the PSA-DME with M^h = 500 
MeV unless otherwise specified. 



V. SELECTED RESULTS 
A. Density-dependent Fock couplings 

In the present section, the non-trivial (isoscalar-) den- 
sity dependence of Fock DME couplings Cf(R;V^^^), 
with g G {/9/9, /9T, yoAp, VpVp, J^p, J J} is briefly ana- 
lyzed. First, we note that C/^''(R; F^N) ^ g through 
N^LO since the two-body spin-orbit interaction is en- 
tirely carried by contact terms (see Eq. IBip . Sec- 
ond, since we are restricting the DME described in 
Sec. mil to the Fock energy contribution, one finds 
^VpVpj.j^.p.NN) ^ Of course, the pAp terms can be 
transformed by partial integration into (Vp)'^ terms, but 
we chose not to do so since the fci?(R)-dependence of the 
couplings results in more complicated expressions. Con- 
versely, we note that it is not possible to convert (Vp)^ 
terms entirely to the pAp form when the couplings de- 
pend on kpCR.) in contrast to usual Skyrme functionals. 

Referring to Figs. [TJIH the non-zero isoscalar and 
isovector couplings are shown including LO, NLO and 
N^LO contributions. The main feature to extract from 
these results is that the non-trivial density-dependence 
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FIG. 1: Density dependence of C/'^(R; V^^) calculated 
through LO, NLO, and N^LO for the isoscalar (upper 
panel) and isovector (lower panel) couplings. A SFR 
mass of Mgfr = 500 MeV was used in the NLO and 
NNLO two-pion exchange potentials. 



is controlled by the longest-range parts of the NN in- 
teraction. Indeed, the density profile of the couplings 
is driven almost entirely by the LO term (one-pion ex- 
change) whereas the NLO and N^LO interactions that 
are built from shorter-range two-pion exchanges provide 
small corrections to the LO density-dependence. One 
might be surprised by the fact that the successive contri- 
butions do not exactly follow the hierarchy LO > NLO 
> N^LO. However, this is somewhat misleading since 
only the Fock contributions from finite-range NN pieces 
are actually shown and cannot be expected to maintain 
such a hierarchy without including N^LO NNN and NN- 
contact contributions. 

The second important result is that all non- 
zero couplings exhibit a substantial isoscalar density- 
dependence^, especially as one goes to small densities. 
Such in-medium dependencies are obviously at variance 
with standard phenomcnological Skyrme parameteriza- 
tions for which only Cf depends on the medium in time- 



^ Note the large y-axis scale used in most of the figures. The rea- 
son that the vertical scales are large becomes clear as one refers 
to typical values of the couplings for standard Skyrme parame- 
terizations; see section [VBI 
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reversal invariant systems. Investigating the impact of 
such non-trivial in-medium dependencies generated by 
pion exchanges is one of the primary long-term goals of 
our project. Of course, many-body correlations gener- 
ated at higher-orders from both short-range-contact and 
long-range-pion physics will provide couplings with ad- 
ditional medium dependencies. As discussed in the next 
section, the impact of such in-medium effects on observ- 
ables and Skyrme phenomenology can only be fully char- 
acterized by performing full-fledged EDF calculations. 



B. Comparison with Skyrme phenomenology 

According to Eq. IHl Fock DME couplings 
Cf(R;V^^) must be complemented with the con- 
tribution C^{A;V^^) whose expressions are given in 
Appendix [B] As already discussed, three additional 
types of contributions to the EDF should be further 
considered. First is the Hartree term given by Eq. [20l 
Although such a contribution solely depends on the 
local part of the density matrix, it does not take the 
form of a local Skyrme-like EDF (Eq. [28]) when treated 
exactly, as advocated here. Second are the Hartree-Fock 
contributions from both the long-range (V^^) and 
short-range (Vj^) parts of the chiral NNN at N^LO. 
The application of the DME to the V^^ contributions, 
which eventually leads to additional density-dependent 
contributions (R; V^^), will be discussed in a separate 
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paper. Last but not least, contributions beyond HF, 
which can hopefuUy be recast into a quasi-local EDF 
form, will add an additional in-medium renormalization 
to the couplings. 

Keeping in mind the above warnings as to what a com- 
plete non-empirical EDF should contain, we now per- 
form a primitive comparison with Skyrme phenomenol- 
ogy with the goal of providing a zeroth-order assessment 
of the non-trivial in-medium dependence of the coupling 
functions. To do so, we compensate for all missing pieces 
beyond Fock DME couplings Cf(R; VjJ^'^) by providing 
an "uncertainty band" generated by imposing natural- 
ness requirements for the coupling constants of the asso- 
ciated bilinear term in the EDF. In the current context, 
naturalness means that the (dimensionless) coupling con- 
stants are of order unity after appropriate combinations 
of the strong interaction scales and are extracted 
from the energy density^. We refer to Ref. [131 for de- 
tails on applying the naive dimensional analysis (NDA) of 
Manohar and Georgi [s^] to Skyrme EDFs. The resulting 
"naturalness band" should only be viewed as an order- 
of-magnitude estimate of the missing higher-order pieces, 
which in any event will carry additional non-trivial den- 
sity dependencies in a fully microscopic approach. Still, 
the fact that most phcnomcnological Skyrme parametcr- 
izations conform to these naturalness bounds [47t pro- 
vides some justification for such a primitive procedure. 

Given that C/^''(R; 1/^^) = through N^LO, we 
do not provide a graphical comparison to Skyrme phe- 
nomenology since such couplings arc entirely provided 
by the "natural" contribution (i.e., the NN spin-orbit 
interaction is a contact interaction). Note, however, 
that the NNN force at N^LO has long-range spin-orbit 
contributions that will provide non-trivial (i.e. density- 
dependent) spin-orbit couplings in the EDF. All other 
isoscalar and isovector couplings are shown in Figs. [S]- 
[HI Compared to the previous section, a logarithmic scale 
is used for the horizontal axis that gives visually more 
weight to lower densities. 

To perform the comparison, we employ a representa- 
tive set of modern Skyrme parametcrizations: SkM* 5l[, 
T22 0, T44 0, TZA g^, SLy4Tsoif SLy5-hT 
Except for SkM*, they all result from recent investiga- 
tions that aimed at pinning down (while keeping the 
overall quality of modern parameterizations) values of the 
tensor couplings C/"' . The Skyrme couplings displayed in 
Figs. OS] are essentially all within the uncertainty band of 
the DME-inspircd coupling functions around saturation 
density. Although only qualitative, this is a very signifi- 
cant, i.e. non-obvious, result. Of course, the rather con- 
servative uncertainty-band used is typically larger than 
the difference between various fine-tuned parameteriza- 
tions. A more systematic non-empirical treatment of all 
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FIG. 5: Fock DME coupling C^-^iR;V^^) augmented 
with a "natural" Skyrme-like contribution (sec text) 
and compared to the corresponding coupling from a 

representative set of Skyrme parameterizations. Upper 
(lower) panel: isoscalar (isovector) coupling. 



contributions to the coupling functions is expected to 
narrow down the uncertainty band and perhaps allow 
one to discriminate between various empirical parame- 
terizations. Of particular interest are tensor couplings 
C/'' whose preferred range^ is not settled yet by the phe- 
nomenology (see Fig. [S]). 

As already stressed, the present status of our approach 
is such that the above comparison should be taken as 
qualitative at best. Only minimal quantitative informa- 
tion can be gleaned from comparisons at the level of the 
EDF couplings as (i) only the total energy is observ- 
able and that (ii) self-consistent effects can be signifi- 
cant, giving rather different final results even if one starts 
with parameterizations whose couplings look alike on the 
scale of Figs. [5][51 Eventually, comparisons based on 
fully self-consistent EDF calculations performed with the 
(scmi-)non-cmpirical energy functionals presented here 
will provide useful quantitative information. Such a 
project is currently underway (36j . 

Lastly, Figs. [9]fT2l compare the various couplings ob- 



We use /tt = 93 MeV and A^^ = 700 MeV when generating the 
NDA estimates. 



The definition of C/'' presently used corresponds to half the one 
of Ref. 0. 
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tained using the original NV-DME and the PSA-DME 
used here. It is reassuring that, while the numerical 
values of the couplings can change somewhat depend- 
ing on which variant of the DME is used, the overall 
density profiles are relatively insensitive to this choice. 
At least within the semi-phenomenological approach out- 
lined in the present paper, the similar density profiles for 
the PSA-DME and NV-DME couplings implies that full- 
fiedged EDF calculations should be fairly insensitive to 
this choice since refitting the Skyrme constants largely 
absorbs these (approximately density-independent) dif- 
ferences. 



VI. SUMMARY AND CONCLUSIONS 

It is by now well-established that empirical Skyrme 
functionals (in present form) exhibit critical limitations 
that are often manifested by parametrization-dependent 
predictions away from known data, e.g., see Ref. (gI. [53|. 
One possible remedy is to use many-body perturbation 
theory and knowledge of the underlying two- and three- 
nuclcon interactions to identify missing microscopic long- 
distance physics, together with density matrix expansion 
techniques to incorporate these missing ingredients into 
existing Skyrme machinery (l2| - [l5| . 

In the present work, our primary focus is approxi- 
mating the spatial non-locality of the Fock contribution 
from the chiral N^LO NN interaction, while keeping the 
Hartree contribution exact and postponing the treatment 
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of the NNN part to a separate paper. The generic struc- 
ture of the chiral interactions is such that each DME 
couphng Cf (t = 0, 1 and g G {pp,pT,pAp,J\/p,JJ}) 
decomposes into a cutoff-dependent couphng constant 
C^{A;V™) arising from the zero-range contact in- 
teractions and a cutoff-independent couphng function 

{Il;V^^) of the density arising from the long-range 
pion exchanges. In the short term, this separation (which 
is completely unambiguous at the Hartree-Fock level) be- 
tween long- and short-distance physics motivates adding 
the (R'jV^^) to existing Skyrme functionals, upon 
which the original Skyrme constants can then be refit 
to data. This semi-phenomenological approach is mo- 
tivated by the observation that the EFT contact terms 
can in principle be fixed to any low-energy quantities. 
Therefore, refitting the Skyrme constants can be viewed 
as matching a microscopically-constrained Skyrme-likc 
functional containing explicit pion physics (albeit in a 
zcroth-order HF approximation at this point) to finite- 
density observables. 

Restricting ourselves to time-reversal invariant sys- 
tems, analytical expressions of the Fock DME couplings 
Cf(R;T/^'^) have been derived as a function of the 
isoscalar density. The rather lengthy analytic expres- 
sions can be downloaded through a companion Mathe- 
matica notebook, and are also collected in Appendix [G] 
The novel density-dependencies of the couplings are con- 
trolled by the longest-ranged parts of the NN interaction, 
which implies that these microscopic constraints are com- 
ing from the best-understood parts of the underlying nu- 
clear interactions. The dependence on the isoscalar den- 
sity is significant for all couplings over the density interval 
of interest, which is obviously at variance with standard 
phenomenological Skyrme paramcterizations whose only 
density-dependent couplings are Cf' {t = 0,1). In the 
longer term, investigating the impact of such non-trivial 
in-medium dependencies generated by pion exchanges is 
one of the central goals of the present project [36l |. 

The rich spin and isospin dependence of the chiral EFT 
one- and two-pion-exchange interactions gives us hope 
that their inclusion will provide valuable microscopic con- 
straints on the poorly-understood isovector properties of 
the EDF. We do not expect dramatic changes for bulk nu- 
clear properties due to the tendency of pious to average 
out in spin and isospin sums, but we do expect inter- 
esting consequences for single-particle properties (which 
phenomenology tells us are sensitive probes of the tensor 
force) and systematics along long isotope chains (which 
should be sensitive to the isovector physics coming from 
pion-exchange interactions). Another potentially signifi- 
cant advantage of the current approach is that two very 
different microscopic origins of spin-orbit properties (i.e., 
short-range NN and long-range NNN spin-orbit interac- 
tions) are treated on equal footing. This is in contrast 
to empirical Skyrme and Gogny functionals, where the 



zero-range spin-orbit interaction has no obvious connec- 
tion with the sub-leading (but quantitatively significant) 
NNN sources of spin-orbit splittings. In a forthcoming 
paper, we will extend the DME to include HF contri- 
butions from chiral three- nucleon interactions at N^LO. 
Such an extension will allow one to probe the impact 
of microscopic threc-nucleon forces on the structure of 
medium- and heavy-mass nuclei. 

The EDF obtained as a result of the present paper and 
the subsequent NNN paper only contains non-empirical 
Hartree-Fock contributions, such that further correla- 
tions must be added to produce any reasonable descrip- 
tion of nuclei. As explained above, in the short term 
this will be done empirically by adding the DME cou- 
plings to empirical and Skyrme functional whose cou- 
pling constants arc refitted to data While this is 
an admittedly empirical procedure, it is motivated by 
the well-known observation that a Brueckner G-matrix 
differs from the vacuum NN interaction only at short 
distances. Therefore, one can interpret the refit to data 
as approximating the short-distance part of the G-matrix 
with a zero-range expansion through second order in gra- 
dients. Eventually though, it is the goal of a future work 
to design a generalized DME that is suited to higher or- 
ders in perturbation theory [s^ . 
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Appendix A: Chiral EFT NN potentials 

The NN potential in chiral EFT can be schematically 
written as 



V 



NN 



= v,T + v,T + --- + Krw, (Ai) 



where Vitt and V2ir arc finite-range one- and two-pion ex- 
change interactions dictated by the spontaneously broken 
chiral symmetry of QCD, and V^^ denotes the scale- 
dependent contact terms that encode the effects of inte- 
grated out degrees of freedom (heavier meson exchanges, 
high-momentum two-nucleon states, etc.) on low energy 
physics. As discussed in the text, our primary focus is 
on the finite-range pion-exchange interactions since they 
drive the non-trivial density-dependencies introduced by 
the DME. In the notation of Eq. [HI we list the non-zero 
finite-range contributions thru N^LO: 
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VT{q) = 



3847r2/^ 
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(g) <^ (5.gl - - 1) + g' 



-^Vs{q) = 

q 



3gi 



L*^=f' {q) 



\2Zg\ - lOff^ - 1) 



48.g>i7r^ 
4toJ + 



(A3) 
(A4) 
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(2m^ + g2)A*^-(a) 



(A5) 
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where the NLO and N^LO loop functions are given by 



^(Aisfr — 2m^) — ni ■ 



2m„ 



2g 4m2(A/4+g2) 

1 g(A/sfi. - 2m^) 

— arctan 

2g g^ + 2A/sfi. m^r 



a; = vV + 4774, 



4m2 , 



(A7) 
(A8) 



Note that in the notations of Equations. [9] and [TOl the 
above finite range interactions do not contribute to the Vi 
and Vf interactions where i G C, S,T. For the numerical 
results presented in the text, we have assumed an axial- 
vector coupling of qa = 1.29, a pion-decay constant of 
/tt — 92.4 MeV, and a pion mass of = 138 MeV. 
The ttttNN low energy constants have been extracted 
from both TriV scattering [U and NN phase shift anal- 
yses [5^, although the C3 and C4 couplings have rather 
large uncertainties whose central values are still some- 
what controversial [HBl • In the results presented here, we 
have adopted the values ci = —0.81 GeV~^, C3 = —3.4 
GeV-\ and C4 = 3.4 GcY'^ as in Ref. [48|. 

In Ref. [131, Epelbaum and collaborators advocate us- 
ing a spectral function regulator (SFR) mass in the range 
of Msir = 500 . . . 800 MeV to reduce the unphysically 
strong attraction in the isoscalar central part of the N^LO 
two-pion exchange potential (TPEP). In contrast, En- 
tem and Machleidt work in a scheme that uses dimen- 
sional regularization to regulate the divergent loop inte- 
grals that enter into the calculation of the TPEP, which 
corresponds to taking Mgh = 00 H^. As emphasized 
in [H^l , formally speaking there is no "correct" value since 
varying Mgh- only modifies the short-distance structure 
of the TPEP, and such variations have no effect on ob- 
servables since they can always be absorbed by the cor- 



r 



responding contact interactions. Nevertheless, it is ar- 
gued in Ref. [l^l that a finite SFR mass in the range of 
A/sfr — 500 . . . 800 MeV offers certain practical advan- 
tages by cutting down the unnaturally strong attraction 
at mid- and short-distances that arises at N^LO. For this 
reason, we have used a value of Afgfr = 500 MeV to 
generate the DME couplings shown in the figures in the 
text. 

On a technical note, we mention that we were not able 
to obtain analytic expressions for the NLO couplings for 
finite values of Mgh ■ However, numerical integration of 
the "master formulas" Eqs. I42ll45] showed rather small 
differences from the Mgh — 00 NLO couplings. For the 
N^LO couplings where the SFR mass has a much more 
quantitative effect, we were able to obtain analytic ex- 
pressions provided Mgh is chosen to be an integer multi- 
ple of m^. For all other choices, such as the Msh = 500 
MeV used in the text, the NLO and N^LO DME cou- 
plings were obtained by numerical integration of Eqs. 

US 



Appendix B: Contact term EDF contributions 

Through N^LO, the EFT two-nuclcon contact interac- 
tion takes the form 
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(PiKriP') 



Cs + Ct (Ti ■ <^2 + CiQ 
+Ce (q • cri)(q • CT2) + C7 (k • CTi)(k ■ aa) , 



C2 + (C3 + C4 k2)(ai . ^72) + iCs i (ai + ^2) • (q x k) 



(Bl) 



where the A-dependence of the coupHngs has been sup- 
pressed for brevity. Decomposing in terms of the 
components {Vc, Vs, . . .} and {Vc, V5, . . .} introduced in 
Eq. [9] gives 



Vc 
Vs 
Vt 
Vls 



Cs + Ciq^ 
Ct + Caq" 



(B2) 
(B3) 
(B4) 
(B5) 



densities multipHed by the density- independent coupHng 
constants in Eas. lB9lfB161 This observation motivates the 
semi-empirical approach advocated in the text whereby 
the density-dependent DME Fock couphngs from the 
finite-range pion exchange contributions are added to ex- 
isting Skyrme functionals, whose empirical parameters 
are then refit to data. In this sense, the refitted Skyrme 
constants can be viewed as containing the contributions 
of Eqs. IB9llB16l in addition to higher order contributions 
beyond HF. 

Appendix C: Exchange interaction 



and 
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By evaluating the action of the exchange operators on 
the direct NN interaction in Eq. [SI the form factors of the 
exchange interaction {Vq, Vg, . . .} appearing in Eq. [TU] 
can be expressed as linear combinations of the direct in- 
teraction {Vc,Vs, ■ ■ ■} as shown in Table [??1 



The EDF couplings arising from Eq. IBll treated at the 
Hartrec-Fock level are given by 
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CPAP ^ l^(^C2 + 3Ci + Cj) 
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C^J = -^(2Ci - 2C3 - 2C4 - 4C6 + Ct) (B15) 
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C^-" = -{Ci-C3-2Ce). 
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TABLE II: Rccoupling coefficients between the direct 

and exchange NN interaction. Note that the 
{Vc,Vs, ■ ■ ■} are evaluated at 2k, since exchanging the 
two nucleons corresponds to replacing q — ?> 2k and 
k ^ q/2. 



As discussed in the text, the Hartree-Fock energy den- 
sity resulting from V^^ is precisely of the same form as 
empirical Skyrme functionals, i.e., bilinear products of 



Appendix D: PSA-DME 



J 



In this appendix we provide a streamlined derivation of the PSA-DME compared to the original presentation in 
R-cf. along with some details concerning the simplified version used in the present work. All variants of the 
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PSA-DME start from the formal identity for the scalar/vector-isoscalar/isovector part of the one-body density matrix 



(ri,r2) = e*'-V'L^^^-'''J ^ (/7*(r2cr2r2) <^»(ri(Tiri) (o-2|<Tp|o-i) (rajf^ 



1 + r • 



i=l 

Vi - V2 



^</J*(r20-2T2) <y9j(ri(Tiri) 



X (o-2|o-^|cri) (r2|'r^|Ti) 



(Dl) 



where indices /.i, e {0, 1, 2, 3}, tq and ctq correspond to a two-by-two identity matrix, ti^2,3 = ''"a;,i/,z and cri,2,3 = 
o'x.y.z- k is a yet-to-be-determined momentum scale to be chosen to optimize the truncated expansion in Eq. IDll 
Physically, k represents an averaged relative momentum in the nucleus. Now, assume we have a model local momentum 
distribution given by g(R, k) and define the following quantities 



n„(r,R) 



/fike*"-'"(r-k)".g(R,k) 



/dk.g(R.k) 

(V«-Vi2))p,.(ri,r2) 



ri=r2=R 



r„6,^,(R) = Vii)vfVA..(ri,r2) 

ri=r2=R 

Performing the phase-space averaging of Eq. (|Dip on the model space defined by g(R., k), we obtain 
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p^^(R) + i 



^ 7-a jpafelR) 



^VaV6p^^(R) - ra6,^^(R) 



(D5) 



where the local densities are as defined previously. 

Even without specifying the actual form of the model momentum distribution, it is clear that the PSA-DME of the 
scalar and vector parts are treated on equal footing (i.e., 11^ = 11^). As shown in Ref. [37[, elementary constraints 
derived from spin-polarized infinite matter forbid the use of channel- dependent Il-functions. Unfortunately, the 
channel-independence of the PSA-DME Il-functions is not entirely transparent from the presentation in Ref. [s^] 
because (i) we only considered the time-reversal invariant case such that only Ilg and ^^'^ to be dealt with for the 
scalar part while only Ilf had to be dealt with for the vector part and (ii) an unnecessary asymmetry was introduced 
in the form of an additional angle-average over the direction of r for the scalar part in order to replace the kinetic 
tensor density Tab with the diagonal kinetic density r. The less transparent derivations in Ref. jssj obscure the fact 
that intrinsically, the PSA leads to a channel- independent DME with the same Il-functions for the scalar and the 
vector parts. In view of the s imp lified (and more general) derivation presented in the current paper, it should be 
realized that the claim in Ref. [331 ^bat the PSA-DME postulates different Il-functions for the scalar and vector parts 
is no longer correct. 

As discussed in Ref. [s^, the PSA-DME is well-suited to incorporate the effects of the diffuseness and anisotropy 
of the local momentum distribution at the spatial surface. However, the inclusion of the diffuseness complicates 
calculations since analytical expressions can no longer be obtained, thus introducing the need for fit parameters into 
the formalism if one desires analytical parameterizations for the couplings. Since the diffuseness primarily affects the 
expansion of the scalar part, which is already reasonably accurate in all existing variants of the DME, we neglect it 
here for simplicity. On the other hand, the inclusion of the anisotropy of the local momentum distribution at the 
nuclear surface predominantly affects the DME of the vector part. As shown in Ref. [s^, the anisotropy does not 
modify the functional form of the Ilf function, as it only enters into the definition of the local Fermi momentum, 
see Eqs. 137] and |3H1 This does not complicate the evaluation of the energy, but it does result in significantly more 
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complicated single-particle fields in the self-consistency loop. Therefore, we take a simplified approach (which still 
gives substantial improvements over the original NV-DME for the vector part) by using the phase space of symmetric 
nuclear matter to perform the averaging, i.e. ^(R, k) = 0(fci? — k). As a result, we find 



noikpr) 
Ui{kpr) 



3^4^«l + 0(fc..r, 
kpr 

-^3Joikpr) + ^9^^ « ^ + 
kpr 5 



U2ikFr) = l5]o{kFr) - 36:^4^ " 3cos(fcFr) « -^^^ + Oikprf . 

kpr 5 
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While Ho starts from 1, the other two H— functions start from 0{kpr)'^. Using a weak ordering argument that counts 
kp on the same ground as the number of gradients, Eq. (jPSP can be rearranged as 



1 , , , , 6abHkpr)kj, ' 



where we neglected i Hi J2a jV'afc(R) that turns out to be a third-order correction (one gradient in the density and 
kp in Ili{kpr)). The prefactor A{kpr) is defined as 



A{kpr) 



^ ^2n,(fc.o + n.(fc..) 

kf. Ilo{kpr) 
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Approximating A{kpr) « 1, one recovers Eq. 44 of Ref. [33] as 



(ri,r2) w nop^^(R) 4- illo ^ ja,/ii/(R) + ^ 
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^VaV6Pp^(R) - Tab,^^(R) 
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(Dll) 



Note that within this approximation scheme (i.e., using the weak fci?-ordering and approximating A{kpr) « 1), the 
constraints on the 11— function resulting from requiring gauge invariance of the energy density functional are satisfied 
trivially. 

In order to recover Eas.[?T]and[5^ one first uses the fact that the one-body density matrix is taken to be diagonal 
in isospin space. This implies that (Op^(ri,r2) is non-zero only if = {0,3}. Hence, we identify the correspondence 



Poo = Po(ri,i"2) 

P03 = Pi(ri,i"2) 

PiQ = so,i(ri,r2) 

Pi3 = si,i(ri,i"2) , 



(D12) 



where ie{l,2,3}. Starting with Eq. (jD12|) one can obtain the corresponding relations involving the local densities. 
For instance, jo,oo(R) is the current density jo,a(R), while ja.io is the tensor spin-orbit density Jo,ai(R)- To obtain 
Eq. pip , we perform angle averaging over the orientation of r. Using the identity 



— / de,. (r • A)(r • B) = —A • B . 



(D13) 



and noting that the current density jQ/x(R) vanishes in time-reversal invariant systems, one obtains Eq. pip with 



Ilg. We reiterate that unlike NV-DME, the IT— functions of PSA-DME satisfy the constraint from gauge 



invariance trivially. As given, Eq. ((32)) corresponds to the DME of St(ri,r2) for spherical systems. To obtain that, 
one combines Eq. (jPlip and Eqs. (|D12p to give 



st,u[ R+|,R-| 



iUlikpr) ^r^Jt,^,(R) 
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Further reduction is possible in spherical systems, for which one can write Jt_^^(R) as a sum of pseudoscalar, vector 
and (antisymmetric) traceless tensor parts 



Jt,iMui'R) - - (5^1. Jt^°^ (R) + - ^ e^ixfe Jt_fc (R) + Jtfj^(R), 

k—x 

where the three components read 
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fl.V — X 

z 

4%{R) ^ J,,^.(R)-i5^.jf'(R) - i jW(R). 



(D15) 

(D16) 
(D17) 
(D18) 



Noting that only the vector part survives in spherical syster ns, one can recovers Eg. ([5^ from Eg. (|D14(). 



Appendix E: Single particle fields 

The density-dependence introduced by the DME of all 
couplings appearing in the EDF makes the single-particle 
fields more complex than the usual case where only the 
Cf couplings arc density-dependent. Note that one can 
generate density-independent couplings by setting the 
DME momentum scale kp{r) egual to a constant, while 
setting fcF(r) = recovers the naive Taylor series expan- 
sion. 

For a systematic comparison of the single-particle 
fields that result from the density-dependent and density- 
independent cases, we consider two cases: one where we 
keep the DME momentum scale intact and the other 
where we set kpir) — 0. Also, we consider a spheri- 
cal system where the single-particle eguation of motion 
is given by 



with the single-particle spinor 



(El) 




The structure of h'^ and its various components read 
/i^ = -V • Br{r)\7 + Ur{r) - iWr • V x ct , 



(E2) 



Ur{r) 



5E 

Spr ' 

SE 



(E3) 



where E ^ E[p, r, J] is the Hartree energy plus the DME 
approximation to the Fock energy. Since in this work the 
density-dependence of the couplings is encoded in fci?(r). 



the only systematic difference between the fields in the 
two cases appears in Ur{r). I.e. 



(E4) 



where U^{r) denotes the field that results after setting 
kpir) = while U!^''{r) is due to the explicit density- 
dependence of the couplings 



_ <5^^|A;jr(r)=0 
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SE 5k F 
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Appendix F: Hartree couplings 

As discussed in the text, it is desirable to treat the 
Hartree energy exactly since the DME is known to per- 
form poorly for such contributions. Nevertheless, it is 
possible to apply the DME to reduce the finite-range 
Hartree energy to the form of a local EDF. For complete- 
ness, we provide expressions that can be used to calcu- 
late the Hartree EDF couplings. In the present case, the 
Hartree energy Eg. [50] simplifies since the only non-zero 
finite-range NN contribution arises from the central force 



Vh^IY. /rfRrfrpi(R + r/2)pt(R-r/2)r^(r) 



Following Negele and Vautherin [301, take 
Pt(R + r/2)pi(R-r/2)«p?(R) 
+ Ir^gikpr) Lt(R)VVt(R) - |Vpt(R)|2 



(Fl) 
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where g{x) = S5j3{x)/2x^. Inserting Eq.|Flinto Eg. [FT] 
gives the following additional EDF couplings from the 
Hartree energy 



<5Cr-ir*(q = 0), (F3) 



6C^^' J q'dqrl{q)h{q/kj.) , (F4) 



where 



hiq) J x'^dxg{x)jo{qx) , (F6) 
= -^-f{M'-i)0{l~q). (F7) 

In Ref. [s^], the authors advocate using a pure Tay- 
lor series expansion for the Hartree contribution. This 
amounts to setting g{x) ~ 1, which results in 

= (F8) 
q dq"^ 

Using this value of hiq) in Eq. (|F4p . one obtains the 
corresponding EDF couplings from the Hartree energy. 



Appendix G: DME couplings from Mathematica 



In this section, we give explicit expressions for the EDF 
couplings that result from the application of NV-DME 
and PSA-DME to the Fock energy from NN chiral EFT 
interaction thru N'^LO. The couplings shown are calcu- 
lated for Mgfr = oo since the analytical expressions are 
more compact than for finite values of the SFR mass. For 
a more complete listing of the EDF couplings (including 
the capability to re-calculate the couplings using finite 
values of Msb ), refer to the Mathematica notebook that 
is provided with this submission. The symbol names for 
the couplings follows the simple rule: DME-type 4- "C" 
-|- EDF-term -I- iso-scalar/iso- vector -|-"noSFR". The re- 
maining Mathematica symbols are self-explanatory. 



NV-DME couplings 



NVCpptOnoSFR = 
I 1 

he gA^ (u^ (-210 (110 cl + 43 c3 - 56 c4) + 6 (23 342 cl + 12 869 c3 - 27 325 c4) 

,73 920 f^7ru^ 

(-6578 cl + 7385 c3 - 14 098 c4) + 2 (-2607 cl + 1357 (c3 + c4) ) + 570 (c3 + c4) 4 
3 u (-80 850 cl + 2079 (2 cl - c3) + 572 (cl - c3 - c4) + 190 (c3 + c4) - 

1155 (c3 - 2 c4) (35 + 2 u^) ) ArcTan [u] ) - 6 (-35 (110 cl + 43 c3 - 56 c4) - 
352 (54 cl + 23 c3 - 40 c4) + 297 (14 cl + 11 c3 - 24 c4) u^) Log[l + u^] ) + 

1 

3 gA^ (4 (21 - 498 - 64 + 16 + 12 u (35 + 4 u^) ArcTan [2 u] ) + 

2048 

3 (-7 + 16 (-8 + 9 u^)) Log[l + 4 u^] ) - 



[4 (-11 025 (-1 + 10 gA^ + 71 gA'') + 525 (-839 + 3014 gA^ + 25 009 gA"*) u 

11468 800 f^7r2u8 I ^ ^ 

4200 (-17 - 6 gA^ + 247 gA^) u'' + 140 (- 1 + gA) (1 + gA) (l33 + 851 gA^) + 
1536 (-1 + gA) (1 + gA) (l + 11 gA^) u^) + 525 Log[l + 2 u + -s/^TTu^j ] 

^4 u a/ 1 + (-21 + 334 + 56 + 16 + 2 gA^ (l05 - 518 + 40 + 48 u^) - 
gA'' (-1491 + 9122 + 648 + 112 u^) ) + 3 (7 - 70 gA^ - 497 gA^ - 
64 (-1 + 6 gA^ + 43 gA^) + 48 (-1 + 2 gA^ + 23 gA^) u^) Log[l + 2 u |u + ^j 1 + 

NVCpptlnoSFR = 
I 1 

he gA^ (-u^ (210 (330 cl + 129 c3 + 56 c4) - 6 (70 026 cl + 38 607 c3 + 27 325 c4) + 

,221 760 TTU^ 

(19 734 cl - 7 (3165 c3 + 2014 c4) ) + 2 (7821 cl + 1357 (-3 c3 + c4) ) + 570 (-3 c3 + c4) 
+ 3 u (242 550 cl - 6237 (2 cl - c3) u'' - 572 (3 cl - 3 c3 + c4) + 190 (-3 c3 + c4) + 
1155 (3 c3 + 2 c4) (35 + 2 u^) ) ArcTan[u] ) - 6 (-35 (330 cl + 129 c3 + 56 c4) - 

352 (162 cl + 69 c3 + 40 c4) + 891 (14 cl + 11 c3 + 8 c4) u^) Log[l + u^] ) + 

1 

gA^ (-4 (21 - 498 - 64 u" + 16 + 12 u (35 + 4 u^) ArcTan[2 u] ) + 

2048 

1 

3 (7 + 16 u= (8 - 9 u^) ) Log[l + 4 u^] ) + 



34 406 400 TT^ 

4 (11 025 (1 - 10 gA^ + 41 gA"*) - 525 (839 - 3014 gA^ + 25 951 gA^) - 4200 (l7 + 6 gA^ + 249 gA^) 
u'' + 140 (-133 - 718 gA^ + 2339 gA^) + 1536 (-1 - 10 gA^ + 59 gA'') u^) + 
525 Log[l + 2 u |u + a/ 1 + | ] |4 u 1 + (-21 + 334 + 56 + 16 + gA^ (-861 + 9854 + 
504 - 368 u^) + 2 gA^ (l05 - 518 + 40 u"* + 48 u^) ) - 3 (-7 + 70 gA^ - 287 gA'' - 
64 (1-6 gA^ + 37 gA^) + 48 (l - 2 gA^ + 25 gA^) u^) Log[l + 2 u [u + ^ 1 + 



NVCpTtOnoSFR 
1 



he 



gA^ m (u^ (60 (110 cl + 43 c3 - 56 c4) - 6 (6754 cl + 3713 c3 - 7870 c4) + 



12 672 7ru^° 

(946 cl - 955 c3 + 1916 c4) + 12 (44 cl - 29 (c3 + c4) ) + 54 (c3 + c4) - 

6 u (-5775 (2 cl + c3 - 2 c4) + 99 (2 cl - c3) u'' + 88 (-cl + c3 + c4) + 90 (c3 + c4) u^) 

ArcTan [u] ) +6 (-1100 cl - 430 c3 + 560 c4 - 
99 (54 cl + 23 c3 - 40 c4) + 99 (14 cl + 11 c3 - 24 c4) u^) Log[l + u^] ) + 

35 gA^ (4 (-3 + 72 + 4 - 60 u ArcTan [2 u] ) +3 (l + 18 - 24 u^) Log[l + 4 u^] ) 

+ 

2048 

[-44 100 (-1 + 10 gA^ + 71 gA'') + 4200 (-421 + 1516 gA^ + 12 569 gA^) + 

6 881 280 f^TT^u^" V 

420 (-283 - 30 gA^ + 4233 gA^) + 1680 (-1 + gA) (1 + gA) (l9 + 125 gA^) + 
464 (l + 10 gA^ - 11 gA^) + 105 Log[l + 2 u + A/"lTu^j ] 

4 u ^J 1 (-105 + 1685 + 134 + 8 - 16 + 10 gA^ (l05 - 527 + 22 + 8 - 16 u 

gA^ (7455 - 11 (4205 + 134 + 8 - 16 u^) ) ) + 105 (l - 10 gA^ - 71 gA^ - 

9 (-1 + 6 gA^ + 43 gA^) + 8 (-1 + 2 gA^ + 23 gA^) u^) Log[l + 2 u [u + + 



NVCprtlnoSFR 
1 

he 



38 016 TTu^ 



gA^ m (u^ (60 (330 cl + 129 c3 + 56 c4) - 6 (20 262 cl + 11 139 c3 + 7870 c4) + 



(2838 cl - 2865 c3 - 1916 c4) + 12 (132 cl + 29 (- 3 c3 + c4) ) + 54 (3 c3 - c4) + 
6 u (5775 (6 cl + 3 c3 + 2 c4) + 297 (-2 cl + c3) + 88 (3 cl - 3 c3 + c4) + 

90 (-3 c3 + c4) u^) ArcTan [u] ) + 6 (-10 (330 cl + 129 c3 + 56 c4) - 
99 (162 cl + 69 c3 + 40 c4) + 297 (14 cl + 11 c3 + 8 c4) u^) Log[l + u^] ) + 
35 gA^ (-4 (-3 + 72 + 4 - 60 u ArcTan[2 u] ) +3 (-1 - 18 + 24 u'') Log[l + 4 u^] ) 

6144 f2 

[44 100 (1 - 10 gA^ + 41 gA'') - 4200 (421 - 1516 gA^ + 13 031 gA^) - 

20 643 840 f ^ tt^ u^" ^ 

420 (283 + 30 gA^ + 4583 gA^) + 1680 (-19 - 106 gA^ + 365 gA^) - 
464 (-1 - 10 gA^ + 59 gA^) u^" + 105 Log[l + 2 u |u + YTTu^j ] 

^4 u a/i + (-105 + 1685 + 134 u'' + 8 - 16 + 10 gA^ (l05 - 527 + 22 u'' + 8 - 16 u 
gA'' (-4305 + 49 825 + 1054 - 472 + 944 u^) ) - 105 (-1 + 10 gA^ - 41 gA^ - 

9 (1 - 6 gA^ + 37 gA^) + 8 (l - 2 gA^ + 25 gA'') u^) Log[l + 2 u |u + 1 + 



23 



NVCpAftOnoSFR 
1 



he 



gA^ m (u^ (-60 (110 cl + 43 c3 - 56 c4) + 6 (6754 cl + 3713 c3 - 7870 c4) + 



50 688 7ru^° 

(-946 cl + 955 c3 - 1916 c4) + 12 (-44 cl + 29 (c3 + c4) ) - 54 (c3 + c4) + 

6 u (-5775 (2 cl + c3 - 2 c4) + 99 (2 cl - c3) u'' + 88 (-cl + c3 + c4) + 90 (c3 + c4) u^) 

ArcTan [u] ) - 6 (-1100 cl - 430 c3 + 560 c4 - 
99 (54 cl + 23 c3 - 40 c4) + 99 (14 cl + 11 c3 - 24 c4) u^) Log[l + u^] ) + 
35 gA^ (-4 (-3 + 72 + 4 u" - 60 u ArcTan[2 u] ) +3 (-1 - 18 + 24 u'') Log[l + 4 u^] ) 

8192 f2 

44 100 (-1 + 10 gA^ + 71 gA^) + 4200 (-421 + 1516 gA^ + 12 569 gA^) + 



27 525 120 f ^ TT^ u^" 

420 (-283 - 30 gA^ + 4233 gA^) + 1680 (-1 + gA) (1 + gA) (l9 + 125 gA^) + 
464 (l + 10 gA^ - 11 gA^) + 105 Log[l + 2 u + A/"lTu^j ] 



4 u ^J 1 (-105 + 1685 + 134 + 8 - 16 + 10 gA^ (l05 - 527 + 22 + 8 - 16 + 

gA^ (7455 - 11 (4205 + 134 + 8 - 16 u^) ) ) + 105 (l - 10 gA^ - 71 gA^ - 

9 (-1 + 6 gA^ + 43 gA^) + 8 (-1 + 2 gA^ + 23 gA^) u^) Log[l + 2 u [u + 1 + u 



NVCpAiflblnoSFR 
1 

he 



gA^ m (u^ (-60 (330 cl + 129 c3 + 56 c4) + 6 (20 262 cl + 11 139 c3 + 7870 c4) + 
152 064 f* TTu" 

(-2838 cl + 2865 c3 + 1916 c4) - 12 (132 cl + 29 (-3 c3 + c4) ) + 54 (-3 c3 + c4) - 
6 u (5775 (6 cl + 3 c3 + 2 c4) + 297 (-2 cl + c3) + 88 (3 cl - 3 c3 + c4) 

+ 90 (-3 c3 + c4) u^) ArcTan[u]) - 6 (-10 (330 cl + 129 c3 + 56 c4) - 
99 (162 cl + 69 c3 + 40 c4) + 297 (14 cl + 11 c3 + 8 c4) u'') Log[l + u^] ) + 
35 gA^ (4 (-3 + 72 + 4 - 60 u ArcTan [2 u] ) +3 (l + 18 - 24 u^) Log[l + 4 u^] ) 

24 576 f2 u^" 

44 100 (l - 10 gA^ + 41 gA'') - 4200 (421 - 1516 gA^ + 13 031 gA^) - 



82 575 360 f ^ tt^ u^" 

420 (283 + 30 gA^ + 4583 gA^) + 1680 (-19 - 106 gA^ + 365 gA^) - 

464 (-1 - 10 gA^ + 59 gA^) u^" + 105 Log[l + 2 u |u + YTTu^j ] 

^4 u a/i + (-105 + 1685 + 134 u'' + 8 - 16 + 10 gA^ (l05 - 527 + 22 u'' + 8 - 16 u^) + 
gA'' (-4305 + 49 825 + 1054 - 472 + 944 u^) ) - 105 (-1 + 10 gA^ - 41 gA^ - 

9 (l - 6 gA^ + 37 gA^) + 8 (l - 2 gA^ + 25 gA'') u^) Log[l + 2 u |u + 1 + 



NVCJJtOnoSFR 



he gA^ m (-4 c4 (l + u^) - 6 cl (l + 3 u^) + 3 c3 (l + 7 + 8 u^) + 8 (l + u^) 

3072 TTu^ (l + u^) 

(2 u (-3 cl + 3 (c3 + c4) +2 (3 c3 + c4) u^) ArcTan[u] - c4 Log[l + u^] ) ) + 

gA^ (4 (5 + 12 u^) + (l + 4 u^)^ Log[l + 4 u^] ) i 

+ 

128 f= (1 + 4 u^) ^ 12 288 (l + u^) 

|-2 a/i + (-1 - 3 - 2 - 4 gA^ (l + 6 + 5 u^) + gA^ (8 + 54 + 52 u^) ) + 

3 Log[l + 2 u + V 1 + j ] ^2 u ( (l + u=)% 2 gA^ (4 + 9 + 5 u^) - gA" (lO + 28 + 17 

gA^ (l + u^)^'^ Log[l + 2 u fu + ^['l7^ 



NVCJJtlnoSFR 



he 



gA^ m (-18 cl + 9 c3 + (-54 cl + 63 c3 + 4 c4) + 4 (18 c3 + c4) 



9216 TTU^ (l + u^) 

16 (9 cl - 9 c3 + 3 c4 + 2 (- 9 c3 + c4) u^) (u + u^) ArcTan [u] + 8 c4 (l + u^) Log[l + u^] ) - 

gA^ (4u2 (5 + 12u2) + (l + 4 u^) ^ Log [l + 4 u^] ) i 

+ 

384 f2 (l + 4 u^)^ 36 864 tt^ (l + u^)^^^ 

2 -s/l (-1 - 3 - 2 - 4 gA^ (l + 6u^ + 5 u^) + gA^ (8 + 30 + 28 u^) ) + 

6 u (- (l + u^)^ - 2 gA^ (4 + 9 + 5 u^) + gA^ (34 + 76 + 41 u^) ) Log[l + 2 u |u + 1 + j 

9gA^ (l+u^)^'' Log[l + 2u |u + Vl + u^ j] 



25 



PSA - DME couplings 



PSACpptOnoSFR 



he 



gA^ (30 (54 cl + 23 c3 - 40 c4) + 3 (1242 cl - 1027 c3 + 1808 c4) 



, 67 200 f^7ru^ 

4 (-1143 cl + 556 (c3 + c4) ) - 744 (c3 + c4) + 6 (l050 (c3 - 2 c4) + 

945 (-2 cl + c3) + 288 (-cl + c3 + c4) u'' - 124 (c3 + c4) u^) ArcTan [u] + 3 (-230 c3 + 
400 c4 + 9 (-60 cl - 12 (14 cl + 11 c3 - 24 c4) + 7 (10 cl - 7 c3 + 8 c4) u^) ) Log[l + u^] ) + 

1 

3 gA^ (4 (15 - 54 + 40 + 120 u ArcTan [2 u] ) - 3 (5 + 24 (3 + u^) ) Log[l + 4 u^] ) + 



5120 

1 



900 (-1 + 6 gA^ + 43 gA^) - 60 (-373 - 258 gA^ + 5095 gA^) + 



1 638 400 TT^ 

320 (31 + 166 gA^ - 197 gA^) - 1104 (-1 + gA) (1+gA) (l + 11 gA^) + 

60 u ^J 1 + (15 - 122 - 80 u'' - 6 gA^ (l5 + 38 + 80 u^) + gA^ (-645 + 1214 + 560 u^) ) 
Log[l + 2 u |u + a/T7u^| ] + 45 (-5 - 24 (2 + u^) + 6 gA^ (5 + 16 - 8 u^) + 



PSACppblnoSFR 
1 



he 



gA^ (215 + 1104 + 264 u^) ) Log[l + 2 u |u + 1 + xi^ j ] 
gA^ (-30 (162 cl + 69 c3 + 40 c4) + 3 (-3726 cl + 3081 c3 + 1808 c4) u'' 



4 (342 9 cl + 556 (-3 c3 + c4) ) + 744 (3 c3 - c4) + 6 (-1050 (3 c3 + 2 c4) - 

2835 (-2 cl + c3) + 288 (3 cl - 3 c3 + c4) + 124 (3 c3 - c4) u') ArcTan[u] + 3 (690 c3 + 
400 c4 + 9 (180 cl + 36 (14 cl + 11 c3 + 8 c4) + 7 ( - 30 cl + 21 c3 + 8 c4) u^) ) Log[l + u^] ) + 

1 

gA^ (-4 (15 - 54 + 40 + 120 u ArcTan[2 u] ) + 3 (5 + 24 (3 + u^) ) Log[l + 4 u^] ) - 



5120 

1 



-900 (1 - 6 gA^ + 37 gA^) + 60 (373 + 258 gA^ + 4777 gA^) - 



4 915 200 7^ 

320 (-31 - 166 gA^ + 533 gA^) - 1104 (-1-10 gA^ + 59 gA^) + 

u ^j (15 (1-6 gA^ + 37 gA^) - 2 (61 + 114 gA^ + 433 gA^) + 80 (-1-6 gA^ + 23 gA^) u^) 



60 ' 



Log[l + 2 u + a/ 1 + j ] - 45 (5 + 24 + + 6 gA^ (-5 + 8 (-2 + u^) ) 

5 gA^ (37 + 24 (lO + u^) ) ) Log[l + 2 u + a/ 1 + j ] 



PSACpttOnoSFR 
1 



he 



3 gA^ m (2 (ll c3 - 24 c4 - (c3 + c4) (l6 + u^) + 14 cl (l + 2 u^) - u (70 cl - 



4480 TTU^ 

35 c3 + 28 (-cl + c3 + c4) + 36 (c3 + c4) u^) ArcTan [u] ) + 
(-28 cl - 22 c3 + 48 c4 + 7 (10 cl - 7 c3 + 8 c4) u^) Log[l + u^] ) + 

9 gA^ ( - 4 + ( 1 + 2 ) Log f 1 + 4 1) i , 

^ ^ ' ^ ^' ' 36 (-l + 2gA2 + 23gA^) + 



256 f^m2u^ 24 576 f^7r2u^ 

48 (-1 + gA) (1+gA) (2 + 13 gA^) + (-4 - 40 gA^ + 44 gA^) + 3 Log[l + 2 u + A/TTu^j ] 

|4 u ^J 1 + (-3 + - 2 u'' + gA^ (6 + 10 - 20 u'') + gA^ (69 - 11 + 22 u"*) ) + 
(3 + 9 - 3 gA^ (2 - 6 + gA^ (23 + 33 u^) ) ) Log[l + 2 u [u + ^J 1 + 



PSACpttlnoSFR = he 



gA^ m (2 (33 c3 + 24 c4 - (3 c3 - c4) (l6 + u^) + 42 cl (l + 2 u^) + 



4480 TTU^ 

u (105 (-2 cl + c3) +28 (3 cl - 3 c3 + c4) + 36 (-3 c3 + c4) u^) ArcTan [u] ) + 
(-6 (14 cl + 11 c3 + 8 c4) + 7 (30 cl - 21 c3 - 8 c4) u^) Log[l + u^] ) - 
3 gA^ (-4 + (l + 2 u^) Log[l + 4 u^] ) i 



73 728 f TT^ 



256 

4 (9 - (24 + u^) - 2 gA^ (9 + 66 + 5 u^) + gA'' (225 + 444 + 59 u^) ) + 12 u 1 + 

(-3 + - 2 + gA^ (6 + 10 - 20 u^) + gA^ (- 75 - 59 + 118 u'') ) Log[l + 2 u |u + ^/TT^j ] 
9 (l + 3 + 5 gA'' (5 + 3 u^) + gA^ (-2 + 6 u^) ) Log[l + 2 u |u + ^[T7^jF^ ] 



PSACpAptOnoSFR 
1 

he 



3 gA^ m (2 (-11 c3 + 24 c4 + (c3 + c4) (l6 + u^) - 14 cl (l + 2 u^) + u 



17 920 f* TTU^ 

(70 cl - 35 c3 + 28 (-cl + c3 + c4) + 36 (c3 + c4) u^) ArcTan [u] ) + 
(28 cl + 22 c3 - 48 c4 - 7 (10 cl - 7 c3 + 8 c4) u^) Log[l + u^] ) - 
9 gA^ ( - 4 + ( 1 + 2 ) Log [ 1 + 4 ] ) 

1024 98 304 



f-36 (-1 + 2 gA^ + 23 gA^) + 



48 (-1 + gA) (1 + gA) (2 + 13 gA^) + (-4-40 gA^ + 44 gA^) + 3 Log[l + 2 u + ] 
|4 u ^[Z7^j^ (-3 + - 2 u'' + gA^ (6 + 10 - 20 u'') + gA^ (69 - 11 + 22 u'') ) + 

(3 + 9 - 3 gA^ (2 - 6 + gA^ (23 + 33 u^) ) ) Log[l + 2 u + ^| 1 + 



PSACpApblnoSFR = 
1 

he 



- gA^ m (2 (-3 (11 c3 + 8 c4) + (3 c3 - c4) (l6 + u^) - 42 cl (l + 2 u^) + 
17 920 TTU^ 

u (-105 c3 + 42 cl (5-2 u^) +4 (3 c3 - c4) (7 + 9 u^) ) ArcTan [u] ) + 
(84 cl + 66 c3 + 48 c4 + 7 (-30 cl + 21 c3 + 8 c4) u^) Log[l + u^] ) + 

3 gA^ (-4 + (l + 2 u^) Log[l + 4 u^] ) i 

+ 

1024 f2m^u^ 294 912 TT^ 

4 (9 - (24 + u^) - 2 gA^ (9 + 66 + 5 u^) + gA'' (225 + 444 + 59 u^) ) + 12 u 1 + 

(-3 + - 2 + gA^ (6 + 10 - 20 u^) + gA^ (-75 - 59 + 118 u'') ) Log[l + 2 u |u + ^iT+u^j ] 
9 (l + 3 + 5 gA^ (5 + 3 u^) + gA^ (-2 + 6 u^) ) Log[l + 2 u |u + ^[T7u^^ ] 



PSACJJtOnoSFR 
1 



he 



- 3 gA^ m (u^ (-28 cl (l + 2 u^) +2 c3 (-11 + 16 + u^) - c4 (l2 + 11 (2 + u^) ) 
8960 77 

2 u (70 cl - 35 c3 + 28 (-cl + c3 + c4) + 12 (3 c3 + c4) u'') ArcTan[u] ) + 
(28 cl + 22 c3 + 12 c4 + 7 (-10 cl + 7 c3 + 4 c4) u^) Log[l + u^] ) + 
9 gA^ (-4 + 8 + Log[l + 4 u^] ) i 

1024 f2m2u^ 49 152 TT^ 

■4 (- 9 + 24 + + gA^ (72 - 273 - 53 u^) + 2 gA^ ( 9 + 66 + 5 u^) ) + 12 u ^| 1 + 



(-3 + - 2 + gA^ (6 + 10 - 20 u^) + gA^ (24 - 35 + 34 u^) ) Log[l + 2 u |u + ^fl 
9 (l + 3 + gA^ (-2 + 6 u^) - gA^ (8 + 15 u^) ) Log[l + 2 u + ^/TTu^j ] 



PSACJJtlnoSFR 
1 



he 



- gA^ m (u^ (-84 cl (l + 2 u^) +6 c3 (-11 + 16 + u^) + c4 (l2 + 11 (2 + u^) ) + 
8960 TTu^ 

2 u (210 cl - 105 c3 - 28 (3 cl - 3 c3 + c4) + 12 (9 c3 - c4) u^) ArcTan[u] ) + 
(84 cl + 66 c3 - 12 c4 - 7 (30 cl - 21 c3 + 4 c4) u^) Log[l + u^] ) - 

3 gA^ (-4 + 8 + Log[l + 4 u^] ) i 

+ 

1024 f2m2u^ 147 456 

■4 (- 9 + 24 + + 2 gA^ (9 + 66 + 5 u'') + gA^ (l80 - 93 + 67 u^) ) + 12 u ^| 1 + 



(-3 + - 2 + gA^ (6 + 10 - 20 u^) + gA^ (60 + 13 + 82 u^) ) Log[l + 2 u |u + ^fl 
9 (l + 3 + gA^ (-2 + 6 u^) - gA^ (20 + 39 u^) ) Log[l + 2 u + a/ 1 + j ] 
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